Roles of RNA scaffolding in nanoscale Gag multimerization and selective protein sorting at HIV membranes.

HIV-1 Gag proteins can multimerize upon the viral genomic RNA or multiple random cellular messenger RNAs to form a virus particle or a virus-like particle, respectively. To date, whether the two types of particles form via the same Gag multimerization process has remained unclarified. Using photoactivated localization microscopy to illuminate Gag organizations and dynamics at the nanoscale, here, we showed that genomic RNA mediates Gag multimerization in a more cluster-centric, cooperative, and spatiotemporally coordinated fashion, with the ability to drive dense Gag clustering dependent on its ability to act as a long-stranded scaffold not easily attainable by cellular messenger RNAs. These differences in Gag multimerization were further shown to affect downstream selective protein sorting into HIV membranes, indicating that the choice of RNA for packaging can modulate viral membrane compositions. These findings should advance the understanding of HIV assembly and further benefit the development of virus-like particle-based therapeutics.

HIV-1 Gag MA domain binds to cardiolipin in a binding mode distinct from virus assemble mediator PI(4,5)P2.

The human immunodeficiency virus type 1 (HIV-1) Gag protein is responsible for facilitating HIV-1 virion assembly and budding. Our study demonstrates that cardiolipin (CL), a component found in the inner mitochondrial membrane, exhibits the highest binding affinity to the N-terminal MA domain of the HIV-1 Gag protein within the lipid group of host cells. To assess this binding interaction, we synthesized short acyl chain derivatives of CL and employed surface plasmon resonance (SPR) analysis to determine the dissociation constants (Kd) for CL and the MA domain. Simultaneously, we examined the Kd of D-myo-phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2 ) derivatives, known to play a crucial role in virion formation. Among all the derivatives, Tetra-C7 -CL exhibited the lowest Kd value (Kd = 30.8 ± 6.9 µM) for MA binding on the CL analog-immobilized sensorchip, indicating a higher affinity. Similarly, the Kd value of Di-C7 -PIP2 (Kd = 36.6 ± 4.7 µM) was the lowest on the PI(4,5)P2 analog-immobilized sensorchip. Thus, Tetra-C7 -CL binds to the MA domain using a distinct binding mode while displaying a comparable binding affinity to Di-C7 -PIP2. This discovery holds significant implications for comprehending the virological importance of CL-MA domain binding, such as its subcellular distribution, including mitochondrial translocation, and involvement in viral particle formation in concert with PI(4,5)P2 . Furthermore, this study has the potential to contribute to the development of drugs in the future.

Conformational transitions of the HIV-1 Gag polyprotein upon multimerization and gRNA binding.

During the HIV-1 assembly process, the Gag polyprotein multimerizes at the producer cell plasma membrane, resulting in the formation of spherical immature virus particles. Gag-genomic RNA (gRNA) interactions play a crucial role in the multimerization process, which is yet to be fully understood. We performed large-scale all-atom molecular dynamics simulations of membrane-bound full-length Gag dimer, hexamer, and 18-mer. The inter-domain dynamic correlation of Gag, quantified by the heterogeneous elastic network model applied to the simulated trajectories, is observed to be altered by implicit gRNA binding, as well as the multimerization state of the Gag. The lateral dynamics of our simulated membrane-bound Gag proteins, with and without gRNA binding, agree with prior experimental data and help to validate our simulation models and methods. The gRNA binding is observed to affect mainly the SP1 domain of the 18-mer and the matrix-capsid linker domain of the hexamer. In the absence of gRNA binding, the independent dynamical motion of the nucleocapsid domain results in a collapsed state of the dimeric Gag. Unlike stable SP1 helices in the six-helix bundle, without IP6 binding, the SP1 domain undergoes a spontaneous helix-to-coil transition in the dimeric Gag. Together, our findings reveal conformational switches of Gag at different stages of the multimerization process and predict that the gRNA binding reinforces an efficient binding surface of Gag for multimerization, and also regulates the dynamic organization of the local membrane region itself.

HIV-1 diverts cortical actin for particle assembly and release.

Enveloped viruses assemble and bud from the host cell membranes. Any role of cortical actin in these processes have often been a source of debate. Here, we assessed if cortical actin was involved in HIV-1 assembly in infected CD4 T lymphocytes. Our results show that preventing actin branching not only increases HIV-1 particle release but also the number of individual HIV-1 Gag assembly clusters at the T cell plasma membrane. Indeed, in infected T lymphocytes and in in vitro quantitative model systems, we show that HIV-1 Gag protein prefers areas deficient in F-actin for assembling. Finally, we found that the host factor Arpin, an inhibitor of Arp2/3 branched actin, is recruited at the membrane of infected T cells and it can associate with the viral Gag protein. Altogether, our data show that, for virus assembly and particle release, HIV-1 prefers low density of cortical actin and may favor local actin debranching by subverting Arpin.

Interaction of TSG101 with the PTAP Motif in Distinct Locations of Gag Determines the Incorporation of HTLV-1 Env into the Retroviral Virion.

Human T-cell tropic virus type 1 (HTLV-1) is known to be mainly transmitted by cell-to-cell contact due to the lower infectivity of the cell-free virion. However, the reasons why cell-free HTLV-1 infection is poor remain unknown. In this study, we found that the retrovirus pseudotyped with HTLV-1 viral envelope glycoprotein (Env) was infectious when human immunodeficiency virus type 1 (HIV-1) was used to produce the virus. We found that the incorporation of HTLV-1 Env into virus-like particles (VLPs) was low when HTLV-1 Gag was used to produce VLPs, whereas VLPs produced using HIV-1 Gag efficiently incorporated HTLV-1 Env. The production of VLPs using Gag chimeras between HTLV-1 and HIV-1 Gag and deletion mutants of HIV-1 Gag showed that the p6 domain of HIV-1 Gag was responsible for the efficient incorporation of HTLV-1 Env into the VLPs. Further mutagenic analyses of the p6 domain of HIV-1 Gag revealed that the PTAP motif in the p6 domain of HIV-1 Gag facilitates the incorporation of HTLV-1 Env into VLPs. Since the PTAP motif is known to interact with tumor susceptibility gene 101 (TSG101) during the budding process, we evaluated the effect of TSG101 knockdown on the incorporation of HTLV-1 Env into VLPs. We found that TSG101 knockdown suppressed the incorporation of HTLV-1 Env into VLPs and decreased the infectivity of cell-free HIV-1 pseudotyped with HTLV-1 Env. Our results suggest that the interaction of TSG101 with the PTAP motif of the retroviral L domain is involved not only in the budding process but also in the efficient incorporation of HTLV-1 Env into the cell-free virus.

HIV-1 Gag targeting to the plasma membrane reorganizes sphingomyelin-rich and cholesterol-rich lipid domains.

Although the human immunodeficiency virus type 1 lipid envelope has been reported to be enriched with host cell sphingomyelin and cholesterol, the molecular mechanism of the enrichment is not well understood. Viral Gag protein plays a central role in virus budding. Here, we report the interaction between Gag and host cell lipids using different quantitative and super-resolution microscopy techniques in combination with specific probes that bind endogenous sphingomyelin and cholesterol. Our results indicate that Gag in the inner leaflet of the plasma membrane colocalizes with the outer leaflet sphingomyelin-rich domains and cholesterol-rich domains, enlarges sphingomyelin-rich domains, and strongly restricts the mobility of sphingomyelin-rich domains. Moreover, Gag multimerization induces sphingomyelin-rich and cholesterol-rich lipid domains to be in close proximity in a curvature-dependent manner. Our study suggests that Gag binds, coalesces, and reorganizes pre-existing lipid domains during assembly.

Structure of the HIV immature lattice allows for essential lattice remodeling within budded virions.

For HIV virions to become infectious, the immature lattice of Gag polyproteins attached to the virion membrane must be cleaved. Cleavage cannot initiate without the protease formed by the homo-dimerization of domains linked to Gag. However, only 5% of the Gag polyproteins, termed Gag-Pol, carry this protease domain, and they are embedded within the structured lattice. The mechanism of Gag-Pol dimerization is unknown. Here, we use spatial stochastic computer simulations of the immature Gag lattice as derived from experimental structures, showing that dynamics of the lattice on the membrane is unavoidable due to the missing 1/3 of the spherical protein coat. These dynamics allow for Gag-Pol molecules carrying the protease domains to detach and reattach at new places within the lattice. Surprisingly, dimerization timescales of minutes or less are achievable for realistic binding energies and rates despite retaining most of the large-scale lattice structure. We derive a formula allowing extrapolation of timescales as a function of interaction free energy and binding rate, thus predicting how additional stabilization of the lattice would impact dimerization times. We further show that during assembly, dimerization of Gag-Pol is highly likely and therefore must be actively suppressed to prevent early activation. By direct comparison to recent biochemical measurements within budded virions, we find that only moderately stable hexamer contacts (-12kBT<∆G<-8kBT) retain both the dynamics and lattice structures that are consistent with experiment. These dynamics are likely essential for proper maturation, and our models quantify and predict lattice dynamics and protease dimerization timescales that define a key step in understanding formation of infectious viruses.

The Rous sarcoma virus Gag Polyprotein Forms Biomolecular Condensates Driven by Intrinsically-disordered Regions.

Biomolecular condensates (BMCs) play important roles incellular structures includingtranscription factories, splicing speckles, and nucleoli. BMCs bring together proteins and other macromolecules, selectively concentrating them so that specific reactions can occur without interference from the surrounding environment. BMCs are often made up of proteins that contain intrinsically disordered regions (IDRs), form phase-separated spherical puncta, form liquid-like droplets that undergo fusion and fission, contain molecules that are mobile, and are disrupted with phase-dissolving drugs such as 1,6-hexanediol. In addition to cellular proteins, many viruses, including influenza A, SARS-CoV-2, and human immunodeficiency virus type 1 (HIV-1) encode proteins that undergo phase separation and rely on BMC formation for replication. In prior studies of the retrovirus Rous sarcoma virus (RSV), we observed that the Gag protein forms discrete spherical puncta in the nucleus, cytoplasm, and at the plasma membrane that co-localize with viral RNA and host factors, raising the possibility that RSV Gag forms BMCs that participate in the intracellular phase of the virion assembly pathway. In our current studies, we found that Gag contains IDRs in the N-terminal (MAp2p10) and C-terminal (NC) regions of the protein and fulfills many criteria of BMCs. Although the role of BMC formation in RSV assembly requires further study, our results suggest the biophysical properties of condensates are required for the formation of Gag complexes in the nucleus and the cohesion of these complexes as they traffic through the nuclear pore, into the cytoplasm, and to the plasma membrane, where the final assembly and release of virus particles occurs.

Dynamics of upstream ESCRT organization at the HIV-1 budding site.

In the late stages of the HIV-1 life cycle, membrane localization and self-assembly of Gag polyproteins induce membrane deformation and budding. Release of the virion requires direct interaction between immature Gag lattice and upstream ESCRT machinery at the viral budding site, followed by assembly of downstream ESCRT-III factors, culminating in membrane scission. However, molecular details of upstream ESCRT assembly dynamics at the viral budding site remain unclear. In this work, using coarse-grained (CG) molecular dynamics (MD) simulations, we investigated the interactions between Gag, ESCRT-I, ESCRT-II, and membrane to delineate the dynamical mechanisms by which upstream ESCRTs assemble templated by late-stage immature Gag lattice. We first systematically derived "bottom-up" CG molecular models and interactions of upstream ESCRT proteins from experimental structural data and extensive all-atom MD simulations. Using these molecular models, we performed CG MD simulations of ESCRT-I oligomerization and ESCRT-I/II supercomplex formation at the neck of the budding virion. Our simulations demonstrate that ESCRT-I can effectively oligomerize to higher-order complexes templated by the immature Gag lattice both in the absence of ESCRT-II and when multiple copies of ESCRT-II are localized at the bud neck. The ESCRT-I/II supercomplexes formed in our simulations exhibit predominantly columnar structures, which has important implications for the nucleation pathway of downstream ESCRT-III polymers. Importantly, ESCRT-I/II supercomplexes bound to Gag initiate membrane neck constriction by pulling the inner edge of the bud neck closer to the ESCRT-I headpiece ring. Our findings serve to elucidate a network of interactions between upstream ESCRT machinery, immature Gag lattice, and membrane neck that regulate protein assembly dynamics at the HIV-1 budding site.

Inositol Hexakisphosphate (IP6) Accelerates Immature HIV-1 Gag Protein Assembly toward Kinetically Trapped Morphologies.

During the late stages of the HIV-1 lifecycle, immature virions are produced by the concerted activity of Gag polyproteins, primarily mediated by the capsid (CA) and spacer peptide 1 (SP1) domains, which assemble into a spherical lattice, package viral genomic RNA, and deform the plasma membrane. Recently, inositol hexakisphosphate (IP6) has been identified as an essential assembly cofactor that efficiently produces both immature virions in vivo and immature virus-like particles in vitro. To date, however, several distinct mechanistic roles for IP6 have been proposed on the basis of independent functional, structural, and kinetic studies. In this work, we investigate the molecular influence of IP6 on the structural outcomes and dynamics of CA/SP1 assembly using coarse-grained (CG) molecular dynamics (MD) simulations and free energy calculations. Here, we derive a bottom-up, low-resolution, and implicit-solvent CG model of CA/SP1 and IP6, and simulate their assembly under conditions that emulate both in vitro and in vivo systems. Our analysis identifies IP6 as an assembly accelerant that promotes curvature generation and fissure-like defects throughout the lattice. Our findings suggest that IP6 induces kinetically trapped immature morphologies, which may be physiologically important for later stages of viral morphogenesis and potentially useful for virus-like particle technologies.

Dual display hemagglutinin 1 and 5 on the surface of enveloped virus-like particles in silkworm expression system.

Rous sarcoma virus-like particles (RSV-LPs) displaying hemagglutinins of H1N1 (A/New Caledonia/20/99) (H1) and H5N1 (A/Vietnam/1194/2004) (H5) of the influenza A virus were produced. The H1 has its transmembrane domain, but the H5 was fused with the transmembrane domain of glycoprotein 64 (BmGP64) from Bombyx mori nucleopolyhedrovirus (BmNPV). H1 and RSV Gag protein were coexpressed in the hemolymph of silkworm larvae, copurified, and confirmed RSV-LP displaying H1 (VLP/H1). Similarly, the RSV-LP displaying H5 (VLP/H5) production was also achieved. Using fetuin agarose column chromatography, RSV Gag protein-coexpressed H1 and H5 in silkworms were copurified from the hemolymph. By immuno-TEM, H1 and H5 were observed on the surface of an RSV-LP, indicating the formation of bivalent RSV-LP displaying two HAs (VLP/BivHA) in the hemolymph of silkworm larvae. VLP/H1 induced the hemagglutination of red blood cells (RBCs) of chicken and rabbit but not sheep, while VLP/H5 induced the hemagglutination of RBCs of chicken and sheep but not rabbit. Additionally, VLP/BivHA allowed the hemagglutination of RBCs of all three animals. Silkworm larvae can produce RSV-LPs displaying two HAs and is a promising tool to produce the bivalent enveloped VLPs for the vaccine platform.

Molecular Determinants of Human T-cell Leukemia Virus Type 1 Gag Targeting to the Plasma Membrane for Assembly.

Assembly of human T-cell leukemia virus type 1 (HTLV-1) particles is initiated by the trafficking of virally encoded Gag polyproteins to the inner leaflet of the plasma membrane (PM). Gag-PM interactions are mediated by the matrix (MA) domain, which contains a myristoyl group (myr) and a basic patch formed by lysine and arginine residues. For many retroviruses, Gag-PM interactions are mediated by phosphatidylinositol 4,5-bisphosphate [PI(4,5)P2]; however, previous studies suggested that HTLV-1 Gag-PM interactions and therefore virus assembly are less dependent on PI(4,5)P2. We have recently shown that PI(4,5)P2 binds directly to HTLV-1 unmyristoylated MA [myr(-)MA] and that myr(-)MA binding to membranes is significantly enhanced by inclusion of phosphatidylserine (PS) and PI(4,5)P2. Herein, we employed structural, biophysical, biochemical, mutagenesis, and cell-based assays to identify residues involved in MA-membrane interactions. Our data revealed that the lysine-rich motif (Lys47, Lys48, and Lys51) constitutes the primary PI(4,5)P2-binding site. Furthermore, we show that arginine residues 3, 7, 14 and 17 located in the unstructured N-terminus are essential for MA binding to membranes containing PS and/or PI(4,5)P2. Substitution of lysine and arginine residues severely attenuated virus-like particle production, but only the lysine residues could be clearly correlated with reduced PM binding. These results support a mechanism by which HTLV-1 Gag targeting to the PM is mediated by a trio engagement of the myr group, Arg-rich and Lys-rich motifs. These findings advance our understanding of a key step in retroviral particle assembly.

Characterization of Bovine Foamy Virus Gag Late Assembly Domain Motifs and Their Role in Recruiting ESCRT for Budding.

A large number of retroviruses, such as human immunodeficiency virus (HIV) and prototype foamy virus (PFV), recruit the endosomal sorting complex required for transport (ESCRT) through the late domain (L domain) on the Gag structural protein for virus budding. However, little is known about the molecular mechanism of bovine foamy virus (BFV) budding. In the present study, we report that BFV recruits ESCRT for budding through the L domain of Gag. Specifically, knockdown of VPS4 (encoding vacuolar protein sorting 4), ALIX (encoding ALG-2-interacting protein X), and TSG101 (encoding tumor susceptibility 101) indicated that BFV uses ESCRT for budding. Mutational analysis of BFV Gag (BGag) showed that, in contrast to the classical L domain motifs, BGag contains two motifs, P56LPI and Y103GPL, with L domain functions. In addition, the two L domains are necessary for the cytoplasmic localization of BGag, which is important for effective budding. Furthermore, we demonstrated that the functional site of Alix is V498 in the V domain and the functional site of Tsg101 is N69 in the UBC-like domain for BFV budding. Taken together, these results demonstrate that BFV recruits ESCRT for budding through the PLPI and YGPL L domain motifs in BGag.

The Characterization of a Novel Virus Discovered in the Yeast Pichia membranifaciens.

Mycoviruses are widely distributed across fungi, including the yeasts of the Saccharomycotina subphylum. This manuscript reports the first double-stranded RNA (dsRNA) virus isolated from Pichia membranifaciens. This novel virus has been named Pichia membranifaciens virus L-A (PmV-L-A) and is a member of the Totiviridae. PmV-L-A is 4579 bp in length, with RNA secondary structures similar to the packaging, replication, and frameshift signals of totiviruses that infect Saccharomycotina yeasts. PmV-L-A was found to be part of a monophyletic group within the I-A totiviruses, implying a shared ancestry between mycoviruses isolated from the Pichiaceae and Saccharomycetaceae yeasts. Energy-minimized AlphaFold2 molecular models of the PmV-L-A Gag protein revealed structural conservation with the Gag protein of Saccharomyces cerevisiae virus L-A (ScV-L-A). The predicted tertiary structure of the PmV-L-A Pol and other homologs provided a possible mechanism for totivirus RNA replication due to structural similarities with the RNA-dependent RNA polymerases of mammalian dsRNA viruses. Insights into the structure, function, and evolution of totiviruses gained from yeasts are essential because of their emerging role in animal disease and their parallels with mammalian viruses.

Visualization of Retroviral Gag-Genomic RNA Cellular Interactions Leading to Genome Encapsidation and Viral Assembly: An Overview.

Retroviruses must selectively recognize their unspliced RNA genome (gRNA) among abundant cellular and spliced viral RNAs to assemble into newly formed viral particles. Retroviral gRNA packaging is governed by Gag precursors that also orchestrate all the aspects of viral assembly. Retroviral life cycles, and especially the HIV-1 one, have been previously extensively analyzed by several methods, most of them based on molecular biology and biochemistry approaches. Despite these efforts, the spatio-temporal mechanisms leading to gRNA packaging and viral assembly are only partially understood. Nevertheless, in these last decades, progress in novel bioimaging microscopic approaches (as FFS, FRAP, TIRF, and wide-field microscopy) have allowed for the tracking of retroviral Gag and gRNA in living cells, thus providing important insights at high spatial and temporal resolution of the events regulating the late phases of the retroviral life cycle. Here, the implementation of these recent bioimaging tools based on highly performing strategies to label fluorescent macromolecules is described. This report also summarizes recent gains in the current understanding of the mechanisms employed by retroviral Gag polyproteins to regulate molecular mechanisms enabling gRNA packaging and the formation of retroviral particles, highlighting variations and similarities among the different retroviruses.

Co-evolution of drug resistance and broadened substrate recognition in HIV protease variants isolated from an Escherichia coli genetic selection system.

A genetic selection system for activity of HIV protease is described that is based on a synthetic substrate constructed as a modified AraC regulatory protein that when cleaved stimulate l-arabinose metabolism in an Escherichia coli araC strain. Growth stimulation on selective plates was shown to depend on active HIV protease and the scissile bond in the substrate. In addition, the growth of cells correlated well with the established cleavage efficiency of the sites in the viral polyprotein, Gag, when these sites were individually introduced into the synthetic substrate of the selection system. Plasmids encoding protease variants selected based on stimulation of cell growth in the presence of saquinavir or cleavage of a site not cleaved by wild-type protease, were indistinguishable with respect to both phenotypes. Also, both groups of selected plasmids encoded side chain substitutions known from clinical isolates or displayed different side chain substitutions but at identical positions. One highly frequent side chain substitution, E34V, not regarded as a major drug resistance substitution was found in variants obtained under both selective conditions and is suggested to improve protease processing of the synthetic substrate. This substitution is away from the substrate-binding cavity and together with other substitutions in the selected reading frames supports the previous suggestion of a substrate-binding site extended from the active site binding pocket itself.

Low CCR5 expression protects HIV-specific CD4+ T cells of elite controllers from viral entry.

HIV elite controllers maintain a population of CD4 + T cells endowed with high avidity for Gag antigens and potent effector functions. How these HIV-specific cells avoid infection and depletion upon encounter with the virus remains incompletely understood. Ex vivo characterization of single Gag-specific CD4 + T cells reveals an advanced Th1 differentiation pattern in controllers, except for the CCR5 marker, which is downregulated compared to specific cells of treated patients. Accordingly, controller specific CD4 + T cells show decreased susceptibility to CCR5-dependent HIV entry. Two controllers carried biallelic mutations impairing CCR5 surface expression, indicating that in rare cases CCR5 downregulation can have a direct genetic cause. Increased expression of ß-chemokine ligands upon high-avidity antigen/TCR interactions contributes to autocrine CCR5 downregulation in controllers without CCR5 mutations. These findings suggest that genetic and functional regulation of the primary HIV coreceptor CCR5 play a key role in promoting natural HIV control.

Rab3a, a small GTP-binding protein, is required for the stabilization of the murine leukaemia virus Gag protein.

We recently identified a CD63-interacting protein to understand the role of CD63 in virion production of the human immunodeficiency virus type 1, and we have found that Rab3a forms a complex with CD63. In this study, we analysed the effect of Rab3a on virion production of the murine leukaemia virus (MLV), which is another member of the retrovirus family. We found that Rab3a silencing induced lysosomal degradation of the MLV Gag protein, and recovery of the Rab3a expression restored the level of the Gag protein through a complex formation of MLV Gag and Rab3a, indicating that Rab3a is required for MLV Gag protein expression. In contrast, CD63 silencing decreased the infectivity of released virions but had no effect on virion production, indicating that CD63 facilitates the infectivity of released MLV particles. Although Rab3a induced CD63 degradation in uninfected cells, the complex of MLV Gag and Rab3a suppressed the Rab3a-mediated CD63 degradation in MLV-infected cells. Finally, we found that the MLV Gag protein interacts with Rab3a to stabilize its own protein and CD63 that facilitates the infectivity of released MLV particles. Considering the involvement of Rab3a in lysosome trafficking to the plasma membrane, it may also induce cell surface transport of the MLV Gag protein.

Quantitative analysis of the formation of nucleoprotein complexes between HIV-1 Gag protein and genomic RNA using transmission electron microscopy.

In HIV, the polyprotein precursor Gag orchestrates the formation of the viral capsid. In the current view of this viral assembly, Gag forms low-order oligomers that bind to the viral genomic RNA triggering the formation of high-ordered ribonucleoprotein complexes. However, this assembly model was established using biochemical or imaging methods that do not describe the cellular location hosting Gag-gRNA complex nor distinguish gRNA packaging in single particles. Here, we studied the intracellular localization of these complexes by electron microscopy and monitored the distances between the two partners by morphometric analysis of gold beads specifically labeling Gag and gRNA. We found that formation of these viral clusters occurred shortly after the nuclear export of the gRNA. During their transport to the plasma membrane, the distance between Gag and gRNA decreases together with an increase of gRNA packaging. Point mutations in the zinc finger patterns of the nucleocapsid domain of Gag caused an increase in the distance between Gag and gRNA as well as a sharp decrease of gRNA packaged into virions. Finally, we show that removal of stem loop 1 of the 5'-untranslated region does not interfere with gRNA packaging, whereas combined with the removal of stem loop 3 is sufficient to decrease but not abolish Gag-gRNA cluster formation and gRNA packaging. In conclusion, this morphometric analysis of Gag-gRNA cluster formation sheds new light on HIV-1 assembly that can be used to describe at nanoscale resolution other viral assembly steps involving RNA or protein-protein interactions.

Genome-Wide Characterization of Zebrafish Endogenous Retroviruses Reveals Unexpected Diversity in Genetic Organizations and Functional Potentials.

Endogenous retroviruses (ERVs) occupy a substantial fraction of mammalian genomes. However, whether ERVs extensively exist in ancient vertebrates remains unexplored. Here, we performed a genome-wide characterization of ERVs in a zebrafish (Danio rerio) model. Approximately 3,315 ERV-like elements (DrERVs) were identified as Gypsy, Copia, Bel, and class I-III groups. DrERVs accounted for approximately 2.3% of zebrafish genome and were distributed in all 25 chromosomes, with a remarkable bias on chromosome 4. Gypsy and class I are the two most abundant groups with earlier insertion times. The vast majority of the DrERVs have varied structural defects. A total of 509 gag and 71 env genes with coding potentials were detected. The env-coding elements were well-characterized and classified into four subgroups. A ERV-E4.8.43-DanRer element shows high similarity with HERV9NC-int in humans and analogous sequences were detected in species spanning from fish to mammals. RNA-seq data showed that hundreds of DrERVs were expressed in embryos and tissues under physiological conditions, and most of them exhibited stage and tissue specificity. Additionally, 421 DrERVs showed strong responsiveness to virus infection. A unique group of DrERVs with immune-relevant genes, such as fga, ddx41, ftr35, igl1c3, and tbk1, instead of intrinsic viral genes were identified. These DrERVs are regulated by transcriptional factors binding at the long terminal repeats. This study provided a survey of the composition, phylogeny, and potential functions of ERVs in a fish model, which benefits the understanding of the evolutionary history of ERVs from fish to mammals. IMPORTANCE Endogenous retroviruses (ERVs) are relics of past infection that constitute up to 8% of the human genome. Understanding the genetic evolution of the ERV family and the interplay of ERVs and encoded RNAs and proteins with host function has become a new frontier in biology. Fish, as the most primitive vertebrate host for retroviruses, is an indispensable integral part for such investigations. In the present study, we report the genome-wide characterization of ERVs in zebrafish, an attractive model organism of ancient vertebrates from multiple perspectives, including composition, genomic organization, chromosome distribution, classification, phylogeny, insertion time, characterization of gag and env genes, and expression profiles in embryos and tissues. The result helps uncover the evolutionarily conserved and fish-specific ERVs, as well as the immune-relevant ERVs in response to virus infection. This study demonstrates the previously unrecognized abundance, diversification, and extensive activity of ERVs at the early stage of ERV evolution.